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INTRODUCTION 


Collisions  between  birds  and  aircraft  have  become  one  of  the  major  flight 
.safety  problems  of  t:lus  jet  ago  •  The  most  critical  need  for  bird  Impact 
resistance  design  Improvement  to  In  military  aircraft'.  Many  high  speed 
aircraft  In  USAP  Inventory  were  not  designed  to  meet  today's  bird  Impact 
reii u l roment a  because  tlio  threat  In  terms  of  lost  aircraft  and  aircrew.*;  was  not 
fully  appreciated  at:  the  design  stage. 

A  comprehensive  review  of  analytically  and  empirically  developed  predic¬ 
tive  techniques  for  bird  impact  structural  performance  has  not  revealed  a 
simple  approach,  that  Is,  simple  to  use  on  short  notice  or  that  requires  a 
minimum  of  Input  information,  nor  one  that  Is  universally  available  and  unl- 
vei  ally  accepted. 

The  statistical  model  for  evaluating  the  bird strike  threat  to  nLrcraft 
transparencies  Is  based  on  two  components'.  (1)  the  operational  Impact  rate 
for  a  specific  usage  profile;  and  (2)  the  probability  that  a  random  hlrdn trike 
wLLl  occur  with  sufficient  kinetic  energy  to  cause  penetration  into  the  air¬ 
craft  cockpit.  This  model  is  based  on  methods  formulated  by  l)r  John  Hal  pin  In 
evaluating  the  birds  trike  threat  on  the  P-16  canopy  and  Dr  Allan  Be  runs'  model 
in  evaluating  the  T-38  crew  enclosure  system. 


This  statistical  simulation  combines  the  velocity  distribution  of  the 
aircraft,  the  bird  weight  distribution,  and  the  capability  of  the  crew  enclos¬ 
ure  expressed  in  terms  of  kinetic  energy.  The  above  distributions  wore  devel¬ 
oped  for  the  K-15  Present  Fleet,  Rapid  Deployment.  Force  and  Dual  Role  Fighter. 
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can  he  expected  on  the  windshield  and  canopy  over  u  realistic  usage  profile. 


The  approach  presented  In  this  report  Is  quite  simple  to  use,  since 
computer  programs  were  generated  to  reduce  the  tedious  calculations  to  a  mini¬ 
mum.  The  use  of  these  methods  will  allow  an  engineer  to  make  a  fairly  quick 
assessment  In  regard  to  the  adcquncy/1  nadequacy  of  a  crew  enclosure  system 
hi rd  res  1  stance  eapahl 1 1 ty . 
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OPERATIONAL  IMPACT  RATE 


As  an  aircraft  flies  through  a  bird  threat  environment,  usually  0  -  5000 
ft  AGL,  it  sweeps  out  a  volume  which  is  a  function  of  the  crew  enclosure  fron¬ 
tal  projected  area,  time  and  aircraft  velocity  in  the  bird  threat  environment. 
Therefore,  It  can  be  concluded  that: 

VOLUME  »  VELOCITY  X  AREA  X  TIME  (1) 

This  volume  per  unit  time  that  Is  swept  out  is  projected  through  the  bird 
threat  environment  with  a  given  bird  density  p  .  The  product  of  the  volume 
per  unit  time  and  bird  density  is  the  Operational  Impact  Rate  (OIR)  on  the 
crew  enclosure  (Figure  1) 


OIR  »  7  X  A  X  Ts  Xp 


(2) 


where 

(HR  *  operational  impact  rate 
V  *  aircraft  mean  true  velocity 
A  =*  windshield  or  canopy  frontal  projected  area 
Ts  «  time  in  bird  threat  environment 
total  flight  time 

p  »  bird  density  in  environment 

A  ratio  can  be  developed  between  any  two  aircraft  if  one  knows  the  operational 
impact  rate  of  one  of  the  aircraft: 


OIR2  ■  OIR j 


V2  X  A2  X  Ts2  X  p2 
Vj  X  Ai  X  Tsj  X  P! 


(3) 


The  mean  velocities,  subsystem  frontal  areas  and  proportions  of  time  are 
directly  measurable  and  can  be  calculated  with  a  moderate  amount  of  effort. 
However,  bird  density  is  not  directly  measurable,  and  if  one  tries  to  deter¬ 
mine  the  bird  density  analytically,  such  a  density  may  be  biased  by  huge  popu¬ 
lations  of  small  birds  such  as  robins  and  blackbirds,  which  have  frequently 
impacted  aircraft.  Another  approach  would  be  to  assume  that  the  two  aircraft 
fly  in  approximately  the  same  environments  and  therefore  experience  the  same 
bird  densities.  This  would  cancel  out  the  two  densities  and  leave  the 
equation  as: 


°Ir2 


OIR  j  X 


V2  x  A2  X  Ts2 
Vj  X  At  X  TS! 


(A) 


Caution  needs  to  be  used  in  this  type  of  approximation,  (assuming  equal  bird 
densities)  since  evaluation  of  existing  bird  impact  data  for  the  F-4  and  F-lli 
(obtained  from  the  Bird  Aircraft  Strike  Hazard  (BASH)  TEAM  at  Tyndall  AFB) 
revealed  a  significant  difference  in  the  predicted  operational  impact  rates. 


2 


OPERATIONAL  IMPACT  RATE  =  (PROJECTED  AREA)  X  (VELOCITY)  X  (TIME)  X  (BIRD  DENSITY) 
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FIG,  1  VOLUME  SWEPT  OUT  BY  AIRCRAFT  IN  BIRD  ENVIRONMENT 


Using  the  known  OIR  for  the  F-4,  the  P-1  It  predicted  OIR  was  high  by  approxi¬ 
mately  a  factor  of  two.  Reversing  the  simulation  showed  the  P-4  predicted  OIR 
to  be  low  by  approximately  a  factor  of  two.  A  realistic  consideration  of  this 
disparity  is  quite  logical,  since  each  aircraft  does  not  fly  at  the  same  exact 
locations  of  the  world  and*  therefore  does  not  encounter  the  same  bird  den¬ 
sities.  Any  extrapolation  between  any  two  unique  aircraft  that  are  not  based 
at  the  same  locations  will  exhibit  this  phenomenon,  and  some  sort  of  correc¬ 
tion  factor  needs  to  be  used. 


In  the  analysis  of  the  F-15  present  fleet.  Rapid  Deployment  Force  (RPD) 
and  Dual  Role  Fighter  (DRF),  the  bird  densities  were  assumed  constant  since 
the  RPD  and  DRF  are  expected  to  fly  at  the  same  bases  from  which  the  histori¬ 
cal  OIR  was  obtained. 


FOR  EXAMPLE:  Windshield,  Europe 


V  X  Ts  X  A 

DRF  DRF  DRF 

OIR  -  OIR  X  - SSESK - 

DRF  Present  ^  X  Ts  X  A 

Europe  Fleet  Present  Present  Present 

Europe  Fleet  Fleet  Fleet 

Europe 


Table  1  summarizes  the  historical  operational  impact  rates  for  the  present 
fleet  based  on  the  time  frame  of  April  1976  -  1  May  1983.  Table  2  summarizes 
the  flight  hours  flown  by  each  F-15  using  command.  Conus  flight  hours  is  a 
summation  of  TAC,  Systems,  and  Logistic  Commands  usage.  Alaskan  Command  and 
PACAF  were  not  analyzed  since  to  date  their  operational  Impact  rate  on  the 
windshield  and  canopy  is  zero,  and  their  hourly  usage  was  minimal  compared  to 
the  usage  in  Conus  and  Europe  (USAFE).  Finally,  Eqn.  (5)  is  used  with  the 
appropriate  inputs  from  Table  3  to  determine  the  predicted  OIR  for  the  RPD  and 
DRF  for  each  of  the  following: 


(A)  Conus  and  Europe 

(B)  Air  to  Air  and  Air  to  Ground  missions 

(C)  Windshield  and  Canopy 

(D)  Ta  variable  from  0  -  .9 
Example  Application  of  Eqn  (5) 

For  F-15  DRF,  Windshield,  Europe,  A/G  Mission  Assume  DRF  t8  ■  .7 

OIR  445.91  .70  602.5 

DRF  -  21.97/106  HRS  X -  X -  X  -  X  1 

EUROPE/ 10^  HRS  279.24  .2307  602.5 

OIR 

DRF  -  106.45  IMPACTS 

EUROPE/ 10*  HRS 
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TABLE  2  F-15  PRESENT  FLEET  TOTAL  USAGE  (FLIGHT  HOURS) 


F-15  PRESENT  FLEET 


ARE  A  i,j  X  Nt)ll  S  l  EL  0 
AREACANOPY 
^AIR  TO  AIR 
CS  CONUS 
CS  EUROPE 


areaWINDSHIELP 
areaCANOPY 
^AIR  TO  AIR 
VAIR  TO  GROUND 
tB  CONUS  AND  EUROPE 


AREAWINDSHIELD 
areaCAN0PY 
VAIR  TO  AIR 
VAIR  TO  GROUND 
cs  CONUS  AND  EUROPE 


TABLE  3 


=.  602.5  IN2 
-  280,0  IN2 

»  279.24  Knots  (CONUS  &  Europe) 

»  .1198  based  on  historical  data 
•»  .2307  based  on  historical  data 


F-15  RAPID  DEPLOYMENT  FORCE  (RPD) 

-  602.5  IN2 

-  280.0  IN2 

-  279.24  Knots  (CONUS  &  Europe) 
■  427.90  Knots  (CONUS  &  Europe) 
«•  Varies  0  -  .9 


F-15  DUAL  ROLE  FIGHTER  (DRP) 

-  602.5  IN2 

-  474.5  IN2 

■  279.24  Knots  ((>6lfoS  &  Europe) 
**  445.91  Knots  (CONUS  &  Europe) 
—  Varies  0  -  .9 


INPUT  VARIABLES  FOR  EQN  (5)  OIR  EXTRAPOLATION 
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Evaluation  of  Table  1  reveals  a  disparity  for  the  Canopy  and  Windshield 
OLRs  between  Conus  and  Europe.  Since  the  canopy  has  a  lower  projected  frontal 
area  from  the  windshield,  it  would  be  expected  that  the  01R  on  the  canopy  in 
Europe  would  be  lower  than  the  windshield  OIR,  however,  this  is  not  the  case. 
Table  4. 0-4. 5  summarizes  the  bird  impacts  for  the  F-15,  F-4  and  F— 111  by  each 
theatre  Conus/Europe  and  each  component  windshield /canopy.  For  Europe,  from 
Table  4.1,  4.3,  4.5  the  results  are  as  follows: 


F-15 

F-4 

F-lll 

Windshield 

3 

38 

20 

Canopy 

7 

60 

15 

A  possible  explanation  of  this  disparity  is  that  a  pressure  wave  is  created  by 
the  aircraft  radome  in  flight,  with  a  component  in  the  vertical  direct  ton 
strong  enough  to  displace  an  oncoming  bird  overA  the  windshield  arch  and  hit 
the  canopy.  The  pressure  force  would  not  be  strong  enough  to  displace  a  4  lb 
bird,  but  It  may  be  sufficient  to  displace  a  .1  .3  lb  bird.  From  Figure  5 
(page  37)  and  6  (page  51)  bird  weight  distributions,  the  bird  weight  range 
0  -  .3  lbs,  represents  35%  of  the  bird  population ‘  in  Conus  and  58%  in  Europe. 

The  above  theory  holds  a  great  deal  of  logical  sense,  but  it  Is  only  a 
crude  rationalization  of  the  results  and  it  falls  for  the  F— 1 11.  An  in-depth 
analysis  is  necessary  to  verify  this  phenomenon  with  special  consideration  in 
the  area  of  aerodynamic  effects  over  the  radome  of  ea :h  of  the  above  aircraft, 
as  well  as  geometrical  similarity.  For  example,  the  F-15  and  F-4  have  basi¬ 
cally  geometrically  similar  windshields  and  canopies  but  the  F-lll  does  not. 
It  is  felt  that  for  the  F-lll  canopy  the  pressure  wave  effects  are  consider¬ 
ably  different  than  the  F-15  and  F-4,  since  the  F-lll  is  a  side-by-side  seat¬ 
ing  aircraft.  Based  on  the  above  the  theory  would  hold  in  Europe  for  the  F-4 
and  F-15,  but  not  for  the  F-lll. 


F-15  CONUS 


ON  CANOPY  0N  WINDSHIELD 


RIKE  DATE 

BASE 

STRIKE  DATE 

BASE 

78/04 

Langley 

78/09 

Unknown 

80/02 

Langley 

79/01 

Robi ns 

79/08 

Holloman 

79/12 

Holloman 

80/02 

George 

80/07 

Langley 

81/07 

Langley 

81/07 

Langley 

81/11 

Langley 

82/05 

Langley 

82/07 

Langley 

83/03 

Holloman 

83/04 

Luke 

83/05 

Luke 

Time  Frame  =  April  1976-May  1983 
Hours  Flown  »  473,008 
Total  -  2  on  Canopy 

14  on  Windshield 

Data  from  BASH  team  computer  files,  Tyndall  AFB 


TABLE  4.0  BIRD  STRIKES  ON  CANOPY  AND  WINDSHIELD 


F— 15  EUROPE 


ON  CANOPY  ON  WINDSHIELD 


STRIKE  DATE 

BASK 

STRIKE  DATE 

BASE 

78/06 

Bitburg 

78/07 

Bitburg 

79/09 

Alconbury 

81/05 

Camp  New  Amsterdam 

80/04 

Camp  New  Amsterdam 

81/07 

Camp  New  Amsterdam 

80/08 

Camp  New  Amsterdam 

82/08 

Spangdahlera 

82/08 

Bitburg 

81/09 

Camp  New  Amsterdam 

Time  Frame  »  April  1976-May  1983 
Hours  Flown  =*  136,576 
Total  =  7  on  Canopy 

3  on  Windshield 


Data  from  BASH  team  computer  files,  Tyndall  AFB 


TABLE  4.1  BIRD  STRIKES  ON  CANOPY  AND  WINDSHIELD 
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F — 111  CONUS 


ON 

CANOPY 

ON  WINDSHIELD 

STRIKE  DATE 

BASE 

STRIKE  DATA 

BASE 

75/07 

Cannon 

75/04 

Nellis 

79/09 

Nellio 

7  5/04 

Cannon 

75/10 

Mt  Home 

76/06 

Mt  Home 

76/08 

Mt  Home 

78/02 

Cannon 

78/06 

Mt  Home 

78/07 

Unknown 

78/06 

Unknown 

78/09 

Unknown 

78/10 

Unknown 

79/08 

Unknown 

78/10 

Unknown 

79/08 

Unknown 

78/10 

Unknown 

79/08 

Unknown 

79/04 

Cannon 

79/08 

Unknown 

79/07 

Unknown 

79/10 

Unknown 

79/10 

Unknown 

79/10 

Unknown 

79/11 

Unknown 

79/11 

Unknown 

79/12 

Unknown 

82/06 

Mt  Home 

82/09 

Mt  Home 

82/09 

Mt  Home 

Time  Frame  =*  Jan  1975-Nov  1982 
Honrs  Flown  =*  Unknown 
Total  =  14  on  Canopy 

16  on  Windshield 

Data  from  BASH  team  computer  flies,  Tyndall  AFH 

TABLE  4.2  BIRD  STRIKES  ON  CANOPY  AND  WINDSHIELD 
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F-lll  EUROPE 


ON  CANOPY  ON  WINDSHIELD 


:rikk  date 

BASK 

STRIKE  DATE 

BASE 

78/07 

Unknown 

75/11 

Upper  lleyford 

78/09 

Unknown 

77/10 

Upper  lleyford 

79/08 

Unknown 

78/08 

Lakenheath 

79/10 

Unknown 

78/08 

Upper  lleyford 

80/07 

Upper  Hey ford 

78/09 

Lakenheath 

80/09 

Upper  Hey ford 

78/09 

Unknown 

81/10 

Lakenheath 

78/09 

Unknown 

82/04 

Incirlik 

78/10 

Unknown 

82/04 

Inclrllk 

79/06 

Unknown 

82/11 

Incirlik 

80/10 

Lakenheath 

82/09 

Upper  Hey  ford 

80/10 

Upper  Hey  ford 

82/10 

Upper  Hey ford 

80/12 

Upper  Hey ford 

82/10 

Upper  Hey ford 

81/01 

Upper  lleyford 

82/10 

Upper  Hey ford 

81/02 

Upper  lleyford 

82/11 

Upper  lleyford 

81/03 

Incirlik 

81/06 

Lakenheath 

81/09 

Upper  lleyford 

81/10 

Unknown 

82/06 

Lakenheath 

82/07 

Upper  lleyford 

Time  Frame  =  , 

Jan  1975-Nov  1982 

Hours  Flown  =» 

Unknown 

Total  »  15  on 

Canopy 

20  oji 

Windshi eld 

Data  from  BASH  team  computer  files,  Tyndall  AFB 

TABLE  4,3  BIRD  STRIKES  ON  CANOPY  AND  WINDSHIELD 


F-4  CONUS 


ON  CANOPY  ON  WINDSHIELD 


strike  date 

BASK 

STRIKE  DATE 

BASE 

77/01 

Boise  (Goweu  Kid)  ID 

75/08 

Holloman 

78/04 

Lincoln  MAP  Nli 

77/01 

Bergstrom 

78/05 

Unknown 

77/03 

Shaw 

79/03 

Unknown 

77/06 

Slutw 

79/04 

Unknown 

77/07 

Birmingham  Muni 

79/04 

Moody 

78/02 

Shaw 

79/04 

Unknown 

78/04 

Seymour  Johnson 

79/04 

Unknown 

78/05 

Unknown 

79/04 

Unknown 

78/05 

Unknown 

79/05 

Moody 

78/06 

George 

79/06 

Unknown 

78/07 

Duluth  IAP  MN 

79/07 

Unknown 

78/08 

Moody 

79/07 

Unknown 

78/08 

Unknown 

79/09 

Unknown 

78/09 

Unknown 

79/12 

Unknown 

78/10 

Unknown 

80/01 

Unknown 

78/11 

Unknown 

80/10 

Duluth  IAP  MN 

78/ It 

Unknown 

81/09 

MacDlll 

78/11 

Moody 

81/01 

Lincoln  MAP  NE 

78/11 

Homestead 

82/07 

George 

79/03 

George 

82/05 

Shaw 

79/05 

Homestead 

82/06 

Birmingham  Muni 

79/06 

Unknown 

82/06 

Hulman  Field  IN 

79/07 

Unknown 

82/11 

Standi  ford  Field  KY 

79/09 

Luke 

82/10 

George 

79/09 

Unknown 

79/10 

Duluth  IAP  MN 

79/10 

Key  Field  MS 

79/10 

Unknown 

79/11 

Homes  tend 

79/1! 

Shaw 

79/12 

Moody 

80/03 

Homestead 

80/04 

Bergstrom 

80/07 

Duluth  IAP  MN 

80/09 

George 

80/11 

Seymour  Johnson 

80/12 

Seymour  Johnson 

81/05 

VDIJi 

81/06 

MacDlll 

81/07 

Bergstrom 

81/09 

George 

81/10 

Kelly 

82/06 

Seymour  .Johnson 

82/08 

Duluth  IAP  MN 

82/10 

Edwards 

82/11 

Seymour  Johnson 

82/  10 

Seymour  Johnson 

82/11 

Seymour  Johnson 

Time  Frame  -  Jail  1975-Jul  1982 
Hours  Flown  “  Unknown 
Total  =  25  on  Canopy 

48  on  Windshield 


Data  from  BASH  team  computer  files,  Tyndall  AFB 


TABLE  4.4  BIRD  STRIKES  ON  CANOPY  AND  WINDSHIELD 


F-4  EUROPE 


ON  CANOPY 


ON  WINDSHIELD 


R.IKK  DATE 

BASE 

STRIKE  DATE 

BASK 

75/03 

Bitburg 

75/09 

Tor re Jon 

76/04 

Alconbury 

78/02 

Bentwatera 

78/03 

Bitburg 

78/03 

Bentwaters 

78/07 

Aviano 

78/06 

Unknown 

78/08 

Torre.1  on 

78/07 

Hahn 

78/09 

Unknown 

78/07 

Unknown 

78/09 

Unknown 

78/08 

Unknown 

78/10 

Unknown 

78/09 

Hahn 

78/10 

Unknown 

78/09 

Unknown 

78/10 

Unknown 

78/09 

Unknown 

78/12 

Zweibrucken 

79/06 

Unknown 

79/01 

Hahn 

79/06 

Unknown 

79/03 

Unknown 

79/06 

Unknown 

79/06 

Unknown 

79/07 

Unknown 

79/06 

Spangdahlem 

79/10 

Hahn 

79/07 

Unknown 

80/07 

Tor re Jon 

79/09 

Unknown 

80/08 

Torrejon 

80/01 

Unknown 

80/08 

Spangdahlem 

80/02 

Unknown 

80/10 

Incirlik 

80/03 

Unknown 

80/11 

Incirlik 

80/0.6 

Torre.1  on 

81/03 

Aleonhury 

80/07 

Zweibruc.ken 

91/03 

Rams tel n 

80/07 

Zweibruc.ken 

81/05 

Spangdahlem 

80/08 

Zaragoza 

81/06 

Rams Coin 

80/09 

Incirlik 

81/06 

Incirlik 

80/09 

Rams te in 

81/07 

Spangdahlem 

80/09 

Zaragoza 

81/07 

Spangdahlem 

80/09 

Hahn 

81/07 

Hahn 

80/09 

Hahn 

81/09 

Torrejon 

80/09 

Zweibruc.ken 

81/09 

Incirlik 

80/09 

Ramstein 

81/10 

Zweibrucken 

80/10 

Zweibrucken 

82/09 

Rams te In 

80/10 

Incirlik 

82/07 

Zweibrucken 

80/10 

Zaragoza 

82/09 

Zweibrucken 

80/10 

Zwoibruo.se  (mi 

82/07 

Torrejon 

80/11 

Raraateln 

82/08 

Torrejon 

81/03 

Incirlik 

82/10 

Torrejon 

81/03 

Incirlik 

81/03 

Incirlik 

81/03 

Ramatein 

81/04 

Incirlik 

81/05 

Incirlik 

81/05 

Hahn 

8  1/07 

Spangdahlem 

TABLE  4.5  BIRD  STRIKES  ON  CANOPY  AND  WINDSHIELD 
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TAIM.K  A.'i  I'1 -A  liUKOPli  (Coat'd) 


81/08 

81/08 

81/09 

81/10 

82/04 

82/05 

82/05 

82/05 

82/06 

82/06 

82/06 

82/09 

82/07 

82/07 

82/07 

82/08 


Torrejon 

Alcoitbury 

Spangdahtem 

Incl  rllk 

Incirllk 

Zaragoza 

Incirllk 

Spaiigdahlem 

Incirllk 

Incirllk 

Incirllk 

Rarastein 

Zweibrucken 

Torrejon 

Torrejon 

Torrejon 


Time  Frame  -  Jan  1975-Jul  1982 
Hours  Flown  *•  Unknown 
Total  -  60  on  Canopy 

38  on  Windshield 

Data  from  KAt.il  team  computer  files,  Tyndall  AFB 
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I 


Cumulative  Probability  Distributions 


Additional  inputs  into  the  model  include  velocity,  bird  weight  and  crew 
enclosure  strength  cumulative  distributions. 

This  model  assumes  that  the  velocity  and  bird  weight  distributions  for  the 
majority  of  the  data  can  be  modeled  by  a  Weibull  cumulative  probability 
distribution. 

The  functional  form  of  the  cumulative  Weibull  distribution  is  given  by: 

(6) 


F(V)  »  1  -  exp 


Where  : 

F(V)  =  cumulative  veiocity  probability  distribution 

V  *  velocity 

6  =  characteristic  paramter,  represents  the  63  percentile  point  of  the 

population  sample 

a  ^  shape  parameter 

Y  =*  minimum  life  parameter,  in  this  analysis  it  is  used  only  as  i 

correction  factor  for  the  Weibull  distribution 

If  portions  of  the  data  population  do  not  fit  a  Weibull  distribution,  an 
approximation  can  be  used  by  the  following  equation 

F(V)fotal  -  F<V)l  +  F(V)2  +  .  .  .  F(V)n  (7) 

Equation  7  sums  each  cumulative  distribution  to  get  a  total  of  1 .  In  this 
report,  F(V)Totap  represents  a  cumulative  probability  distribution  function 
which  is  a  mixture  of  Weibull  distributions  and  distributions  that  are  defined 
by  sixth  or  first  order  equations,  depending  on  which  data  population  is  being 
analyzed. 

The  mean  values  may  be  obtained  by  taking  the  first  derivative  of  each 
individual  cummulative  distribution,  and  applying  the  following  formula 

B 

VL  =  /V  f  ( V) !  dV  (8) 

A 


Where  : 

V  =  velocitv 

f(V>  =  first  derivative  of  F(V)j  probability  distribution 


and 
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VTotal 


(9) 


*  Vj  +  V2  +  Vj  +  .  .  • 

_  B  CD 

vTotal  *  /Vf<V)!  dV  +  /  Vf (V) 2  dV  +  /Vf(V)3  dV  f  .  ..  .  (10) 

A  B  C 

Before  a  cumulative  distribution  can  be  assigned  to  a  given  set  of  data, 
it  Is  necessary  to  reduce  the  data  In  a  statistical  manner.  This  can  be 
accomplished  by  the  following  eqn. 

RELATIVE  FREQUENCY  *  -jt-  (11) 


where 


f  =  frequency  of  occurrence 
N  =  total  number  of  occurrences 

Eqn.  11  suggests  that  one  can  associate  with  a  given  event  a  number,  say 
p,  that  is  equal  or  approximately  equal  to  that  number  about  which  the  rela¬ 
tive  frequency  seems  to  stabilize.  The  number  p,  associated  with  a  given 
event  is  usually  called  a  probability.  Therefore 

p  »  PROBABILITY  *  RELATIVE  FREQUENCY  -  -X-  (12) 


Once  a  given  set  of  data  is  reduced  to  a  probability  form,  it  is  arranged 
in  increasing  order  to  obtain  a  cumulative  probability  distribution.  To  this 
data,  a  known  cumulative  distribution  function  is  fitted.  This  is  done  since 
the  statistically  reduced  data  sample  contains  some  level  of  inaccuracy.  The 
level  of  inaccuracy  is  a  function  of  the  source  and  the  amount  of  data.  For 
example,  structural  recorder  data  is  more  realistic  and  accurate  in  analyzing 
the  F-15  Air  to  Air  mission  than  obtaining  information  from  a  structural  loads 
report . 

A  theoretical  cumulative  distribution  such  as  the  Weibull,  Gama,  Beta  or 
dozens  of  others,  which  can  be  obtained  from  textbooks  are  exact,  meaning  they 
satisfy  all  axioms  of  probability.  If  a  set  of  data  fits  any  theoretical 
distribution  exactly,  the  error  in  that  data  sample  is  zero.  This  occurrence 
is  quite  remote  due  to  the  inaccuracy  in  data  measuring  equipment.  Therefore, 
the  data  is  fitted  with  a  known  cumulative  distribution,  because  such  a 
distribution  is  more  exact  and  represents  actual  aircraft  usage  more  realisti¬ 
cally  for  a  particular  mission  profile. 

This  analysis  uses  primarily  the  Weibull  cumulative  probability  distribu¬ 
tion  since  this  distribution  was  used  exclusively  by  Dr  Halpin  and  Dr  Berens 
in  their  work  in  evaluating  the  F-16  and  T-38  crew  enclosure  bird  impact 
resistance  capabilities.  However,  one  should  note  that  other  known  cumulative 
probability  distributions  could  have  been  used  and  possibly  provided  a  better 
fit  to  any  one  of  the  data  samples  in  this  report. 
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For  the  present  fleet  air  to  air  mission,,  the  velocity  distribution  was 
derived  from  structural  recorder  data  from  Kadena,  Holloman,  and  Bitburg  for  a 
total  of  347.5  flight  hours.  This  data  was  obtained  from  McDonnell  Douglas 
and  reduced  for  an  altitude  range  of  0  -  5000  ft  AGL.  Table  5.0  summarizes 
the  present  usage  of  the  F-15  fleet. 

Eqn  12  was  then  used  to  get  the  initial  probabilities  in  the  altitude 
range  0  -  5000  ft  (Third  Column,  Table  4.1) 

Where 

f  =  ft  seconds  for  specific  airspeed  range 

N  =  total  seconds  (490,981.8) 

From  this  data,  a  cumulative  probability  was  obtained  (Fourth  Column,  Table 
5.1).  This  data  was  then  plotted  on  linear  graph  paper  and  fitted  with  a 
Weibull  cumulative  probability  distribution. 
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DATA  FROM  KADENA,  HOLLOMAN  AND  BITBURG 
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TABLE  5.1  REDUCED  DATA  FOR  0-5000  FT.,  FROM  TABLE  5 


To  determine  the  initial  parameters  for  the  Weibull  distribution,  eqn  (6), 
a  U  vs  W  transformation  was  used.  This  allowed  for  an  initial  estimate  of  the 
Weibull  parameters,  ot  ,  and  3.  The  transformation  was  developed  as  follows: 


F(V)  =  1  -  exp 


1  -  F(  V)  = 


In  (l-F(V)) 

-In  (l-F(V)) 


a 


a-F(v)rj 


a  In  (V  - y)  -  aln  (3 -y)  “  In  Q-  In  (1-F(V))^J 
But  the  equation  of  a  straight  line  is 


u  =  mw  +  b 


where  m  is  the  slope  of  the  line 

b  is  the  u  intercept  when  W  =<  0 


let 

W  = 

In 

(V  -Y) 

and 

u  = 

In 

j~“ln  (l-F<V)r[ 

Therefore ; 

u  = 

m 

W  ■ 

4-  b 

ll  “ 

a 

w 

-  al.n(  3  -  Y  ) 

(13) 


where  a  is  the  slope 

and  -  (x  In  (  3  -  y  )  is  the  u  intercept  when  W  =  0 
*  Any  Weibull  function  in  a  u  vs  W  domain  is  a  straight  line. 

Table  5.2  lists  the  results  of  the  u  vs  V  transformation  where  the  nominal 
velocities  were  input  for  V,  and  cumulative  probabilities,  from  column  3,  were 
input  for  F(V).  These  values  were  then  plotted  in  the  u  vs  W  domain,  (Fig  2) 
and  a  straight  line  was  drawn  through  the  most  significant  points  that  repre¬ 
sented  individual  probabilities,  from  eqn  12,  of  greater  than  .05.  From  this 
straight  line  the  parameters  were  derived  as  follows: 
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-15  PRESENT  FLEET 
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TABLE  5.2  WE I BULL  U  VS  tf  TRANSFORMATION 


a*  slope  of  line  **  2.81 
from  graph,  u  ■  -15.7002  when  W  ■  0 

Y  a  40  established  prior  to  u  vs  W  transformation 


Then 


-15.7002  »  2.81(0)  -2.81  In  (  3  -40) 


3  -  307 


The  correction  factor  (y)  in  the  u  vs  W  domain  straightens  the  data 
somewhat  to  represent  a  straight  line.  The  value  for  Y  was  determined  by 
trial  and  error,  as  y  was  incremented  from  0  -  140,  and  then  finally  chosen  as 
40,  since  it  provided  the  best  fit  to  the  data  sample. 

The  above  Weibull  distribution  fitted  the  data  sample  exceptionally  well 
from  200  to  720  knots,  but  was  unrealistic  below  200  knots.  Therefore,  the 
range  110  -  200  was  fitted  with  a  sixth  order  equation. 

This  distribution  was  developed  by  approximating  that  range  by  a 
realistic  distribution  and  intersecting  the  Weibull  at  the  most  common  point, 
which  for  this  case  was  200  knots.  Then  from  this  distribution,  a  sixth  order 
least  squares  fit  was  accomplished,  and  the  equation  checked  to  Insure  that 
the  distribution  for  the  specified  range  was  in  fact  a  probability  distribu¬ 
tion.  The  following  equation  was  developed: 

(-5.793168  +  2.49687  X  10'1  V  -  4.296046  X  10-3V2 
+3.774985  X  10“5V3  -1.792476  X  10-V*  for  1 16  <  V  <_  200 

+  4.396387  X  lO-l^V3  -4. 371*59  X  10“13v6 

Figure  3  summarizes  the  F-15  present  fleet  velocity  profile  in  the  altitude 
range  of  0  -  5000  ft  AGL. 
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FIG.  2  U  VS  W  TRANSFORMATION 
-IS  PRESENT  FLEET, AIR  TO  AIR  MISSION 


50.0  200.0  250.0  300.0  350.0  400.0  450.0  500.0  550.0  600.0  650 

AIRCRAFT  VELOCITY  (KTS) 

FIG.  3  F-15  PRESENT  FLEET  VELOCITY  DISTRIBUTION 
AIR  TO  AIR  MISSION  FROM  STRUCTURAL  RECORDER  DATA 


The  mean  velocity  was  obtained  by  applying  eqn  9  and  10.  By  using  a 
numerical  integration  subroutine  on  an  HP-15C,  the  mean  velocity  for  the  F-13 
present  fleet  was  determined  to  be  279.24  knots. 

The  proportion  of  time  (tg)  in  the  range  of  0  ~  5000  ft  AGL  for  CONUS  was 
obtained  from  McDonnell  Douglas  structural  loads  report  (MDC  A5318)  for  all 
missions  other  than  air  to  ground.  The  ts  was  determined  to  be  .1198,  based 
on  mission  mix  as  specified  on  pg  3.2.2. 1  of  the  above  report.  The  reason  for 
not  using  the  Holloman  ts  was  that  it  was  unrealistic  and  it  contradicted  a 
data  sample  from  Langley.  The  tg  for  Holloman  was  approximately  .04  and 
Langley,  .14.  Therefore,  as  a  last  resort,  the  structural  loads  report  was 
used. 


Bitburg  data  was  used  to  determine  the  t9  for  Europe.  From  Bitburg  data 
for  Europe,  the  ts  was  evaluated  to  be  .2307.  This  represented  a  realistic, 
value  since  historically  the  F-15  fleet  has  spent  considerably  more  time  in 
Europe  between  0  -  5000  ft,  and  another  source  of  data  matched  the  above  ts 
within  .01  (1%). 

An  assumption  was  made  that  the  velocity  profile  for  the  present  fleet,  in 
the  range  of  0  -  5000  ft  AGL,  was  the  same  between  CONUS  and  Europe.  This  was 
done  considering  that  the  mean  velocity  would  not  increase  or  decrease  by  more 
than  10  knots.  Such  a  change  in  this  analysis  represents  insignificant 
changes  in  final  results. 


F-15  RAPID  DEPLOYMENT  FORCE  (RPD)  VELOCITY  DISTRIBUTION 


For  the  F-15  RPD,  the  air  to  ground  velocity  profile  was  derived  from 
McDonnell  Douglas  structural  load  report  (MDC  A5318).  This  data  was  used 

since  the  F-15  does  not  fly  an  air  to  ground  mission  at  this  time.  The  above 
report,  although  not  labeled  F-15  RPD,  is  representative  of  the  RPD,  since  the 
RPD  will  be  the  basic  F-15C  aircraft  with  provisions  for  conformal  fuel  tanks 
and  a  requirement  for  an  air  to  ground  mission  (based  on  conversations  with 
McDonnell  Douglas). 

Table  6. 0-6. 3  summarizes  the  Air  to  Ground  missions  for  the  F-15  RPD. 

Similar  methods  were  used  to  reduce  the  data  as  in  the  previous  section; 
however,  some  reorganization  of  data  was  necessary  to  fit  the  data  with  a 
Wiebull  curve.  For  example,  if  a  data  point  was  reduced  to  be  at  443  knots, 

then  this  point  was  assigned  to  a  velocity  range  of  410  -  460  knots  and  for 

Weibull  fitting  purposes  the  nominal  point  between  410  -  460  was  used.  This 
approach  allowed  for  some  inaccuracy  in  the  source  of  data.  To  these  bound¬ 
aries,  the  Weibull  distribution  could  have  been  fitted,  iterating  as  necessary 
to  Insure  most  Important  points  are  intersected  or  very  closely  approached. 
Important  points  for  this  case  are  defined  as  points  for  which  eqn  12  shows 
values  of  .15  or  greater  (this  value  can  be  established  by  using  one's  judg¬ 
ment). 

Figure  4  summarizes  the  F-15  RPD  velocity  profile  for  an  altitude  range  of 
0  -  5000  ft  AGL.  An  assumption  was  made  that  the  same  velocity  profile  is 
flown  between  CONUS  and  Europe.  Further,  it  was  assumed  that  for  0  -  5000  ft 
AGL,  the  F-15  RPD  will  fly  the  same  air  to  air  velocity  profile  in  CONUS  and 
Europe  as  the  present  fleet  (Fig  3).  Same  rationale  as  previous  section. 
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PAYLOAD:  4  AIM-7 F  +  Internal  Gun  +  Aa.no  +  4000  lbs.  Bombs 

FUEL:  Full  Internal  ?E?  +  Full  External  £  600  Gal.  Tank. 
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TABLE  6.2  F-15  RAPID  DEPLOYMENT  FORCE  AIR-TO-GROUND  2A  MISSION 


PAYLOAD:  4  AIM-7F  +  Internal  Gun  +  Anno  +  6000  LbSo  Bombs  (MK-82) 


AIRCRAFT  VELOCITY  (KTS) 

FIG.  4  F-15  RPD  VELOCITY  DISTRIBUTION 
AIR  TO  GROUND  MISSION  FROM  STRUCTURAL  LOAD  REPORT 


F-15  DUAL  ROLE  FIGHTER  VELOCITY  DISTRIBUTION 


The  F-15  Dual  Role  Fighter  (DRF),  air  to  ground  (0  -  500  ft  AGL)  velocity 
profile  was  derived  from  McDonnell  Douglas  Structural  Design  Criteria, 
Memorandum  199-DRF-169.  This  report  was  used,  since  it  was  the  best  source  of 
this  type  of  information  for  the  DRF. 


Mission  tnix  was  obtained  from  the  same  report  and  is  summarized  as 
follows : 


F-15  DRF  MISSION  MIX 


AIR-TO-AIR 

20* 

A/Al 

2.5% 

A/A2 

2.5% 

A/A3 

15.0% 

AIR-TO-GROUND 

75* 

A/Gl 

10.0% 

A/G2 

10.0% 

A/G3 

55.0% 

NON-TACT ICAL 

5* 

N/Tl 

2.5% 

N/T2 

2.5% 

NOTE  (A)  Based  on  a  4000  hour  aircraft  lifetime 

(B)  All  percentages  are  percent  of  total  missions 

Table  7.0  -  7.2  lists  each  air  to  ground  mission.  In  data  reduction,  only 
flight  in  the  altitude  range  of  0  -  5000  ft  is  considered.  A  ratio  is  used 
for  climbs  and  descents  up  to  5000  ft.  Same  methods  were  used  as  in  the  pre¬ 
vious  sections  to  fit  the  data  with  a  Weibull  and  sixth  order  distribution. 

Figure  5  summarizes  the  F-15  DRF  velocity  profile  for  an  altitude  range  of 
0  -  5000  ft  AGL.  Similarly  as  for  the  F-15  RPD,  it  was  assumed  that  the  DRF 
will  fly  the  same  velocity  profile  (0  -  5000  ft)  between  CONUS  and  Europe. 
The  DRF  air  to  air  velocity  profile  (0  —  5000  ft)  was  assumed  to  be  equivalent 
to  the  F-15  present  fleet  between  CONUS  and  Europe.  Same  rationale  as  pre¬ 
vious  sections. 
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PAYLOAD:  2  AIM-7F  +  Internal  Gun  +  Ammo  +  8000  lbs.  Bombs  +  4  AIM-9  +  LANTIRN  Pods 
FUEL:  Full  Internal  +  Full  CFT 

ALTITUDE  WEIGHT  FUEL  DISTANCE  TIME 
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PAYLOAD:  2  AIM-9  +  SUU-20  +  LANTIRN  Pods 

FUEL:  Full  Internal  +  Full  OFT 
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VELOCITY  (KTS) 

FIG.  5  F-15  DRF  VELOCITY  DISTRIBUTION 
RIR  TO  GROUND  MISSION  FROM  STRUCTURAL  LORO  REPORT 


BIRO  WEIGHT  DISTRIBUT IONS 


Bird  weight  distributions  for  CONUS  and  Europe  were  derived  from  F-15,  F-4 
and  F— 1 II,  (as  well  as  other  USAF  inventory  aircraft,  Table  8.0)  bird  impacts 
which  were  identified  by  species  and  weight.  The  above  data  was  obtained  from 
Lt  Will  (AV  970-6243),  BASH  Team,  Tyndall  AFB. 

The  bird  impact  data  is  summarized  in  Tables  8.0  -  8.6.  Table  8.0  lists  a 
composite  of  several  aircraft  bird  impacts  which  were  identified  by  the 
Smithsonium  Institute.  Some  identified  bird  weights  are  exact  and  some  are 
averaged.  F-15,  F-4,  and  F-lll  impacts  were  extracted  from  the  log  book  data 
and  listed  in  Table  8.1  -  8.6  with  remaining  computer  file  bird  impacts.  All 
bird  impacts  are  impacts  that  have  been  reported  anywhere  on  the  aircraft.  A 
bird  impacting  a  wing  or  engine  could  have  just  as  well  hit  the  windshield  or 
canopy.  The  data  from  computer  files  was  identified  by  bird  type,  then 
Reference  3-5  was  used  to  associate  a  bird  type  with  a  bird  weight.  Reference 
3  was  considered  to  be  the  most  exact  source  of  information,  then  reference  4 
and  finally  Reference  5  as  a  last  resort. 

Based  on  Tables  8.0  -  8,6  for  CONUS  and  Europe,  the  data  was  reduced  using 
eqn  12  and  then  fitted  with  a  Weibull  distribution  or  a  mixture  of  distribu¬ 
tions.  Figures  6  and  7  summarize  the  results  with  appropriate  equations 
defining  each  distribution,  as  well  as  the  mean  bird  weights. 
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COMPOSITE  OF  SEVERAL  AIRCRAFT  CONUS 
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TABLE  8.0  LIST  OF  IDENTIFIED  BIRDS  INVOLVED  IN  BIRD  STRIKE 


COMPOSITE  OF  SEVERAL  AIRCRAFT  CONUS  (Continued) 
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TABLE  8.0  LIST  OF  IDENTIFIED  BIRDS  INVOLVED  IN  BIRD  STRIKE 
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ims  8.6  LIST  OF  IDENTIFIED  BIRDS  INVOLVED  IN  BIRD  STRIKE 


Ain i auacmd  3-Auynnwna 


BIRD  WEIGHT  (LB. ) 

FIG.  6  CONUS  BIRD  WEIGHT  DISTRIBUTION 
BIRD  WEIGHT  DfiTR  FROM  BfiSH  TERM  COMPUTER  PILES 


umauaoad  ^amutmio 


O  5.0  5.0  7.0  8.0 

BIRD  HEIGHT  (LBS. ) 

EUROPE  BIRD  WEIGHT  DISTRIBUTI 
DOTH  FROM  BASH  TERM  COMPUTER 


COM  PUNK  NT  STRENGTH  DISTRIBUTIONS 


The  canopy  and  windshield  strength  distributions  are  the  percent  area  of 
the  component  which  would  be  damaged,  expressed  as  a  function  of  impact  kine¬ 
tic  energy.  Since  the  majority  of  USAF  bird  impact  testing  Is  done  with  an 
equivalent  of  a  4  lb  bird,  the  strength  distributions  have  been  derived  based 
on  the  kinetic  energy  developed  by  this  bird  weight. 

For  example: 

K.E.  -  1/2  mV2  (14) 

so,  for  the  impact  by  a  4  lb  bird  at  410  Knots: 

m  “ _ 4  ■  .1242  LBS 

32.2  FT/SEC2 

410  Knots  X  1.688  =  692.0800  ft/sec 

K.E.  -  1/2  (.1242)  (692. 0800) 2 

K.E.  -  29,744.3305  FT-LBS 

Based  on  Figure  16  for  the  above  kinetic  energy,  100%  of  the  present  F-15 
windshield  is  critical,  meaning  that  a  penetration  would  occur  anywhere  on  the 
windshLeld  by  a  4  lb  bird.  At  a  velocity  of  330  knots,  an  impact  by  a  4  lb 
bird  (KE  *»  19,269.23  FT-LBS)  makes  the  windshield  44%  critical  (44%  chance  of 
penetration  occurring).  Figure  8  shows  the  F-15  single  and  two  seat  model 
aircraft  windshield  and  canopy  profiles.  Figures  9  and  10  define  the  present 
component  capability  In  terms  of  penet rat lon/no  penetration. 

The  present  capability  for  the  windshield  and  canopy  was  obtained  from 
past  F-15  bird  impact  tests  (McDonnell  Douglas  Report,  MDC  A4888).  Since  two 
bird  impact  [joints  tested  were  not  the  most  critical  locations  on  the  wind¬ 
shield,  T-38  data  was  used  to  develop  the  present  capability  as  shown  in 
figures  9  and  10.  The  T-38  data  was  extracted  from  AFWAL  report  #TR-80-3l  32, 
PART  l,  and  it  was  used  since  the  T-38  had  an  exhaustive  bird  impact  test 
program,  as  well  as  the  T-38  windshield  is  geometrically  similar  and  made  out 
of  the  same  material  (monolithic  stretched  acrylic)  as  the  F-15  windshield. 
Percentage  of  area  was  determined  by  using  the  T-38  percentage  distribution 
(F-15  capability).  Figure  15  shows  the  T-38  penetration  values  by  a  4  IS  bird 
as  obtained  from  the  bird  impact  test  program.  Using  similar  ratios,  the  F-15 
capability  was  obtained  and  is  summarized  in  Figures  9  and  10. 

Example  calculation: 

From  F-15  birdstclke  test  on  windshield,  no  penetration  occurred  at 
393  kts  (severe  cracking);  at  431  kts,  penetration  occurred  (center  impact). 
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The r« fore,  431-393  -  38  KTS 


_38  -  19 

2  393  +  19  -  412  sf  410  KTS 

410  kts  represents  the  predicted  penetration  value  since  severe 
cracking  occured  at  393  KTS. 

* 

=  390  KTS  represents  no  penetration  by  a  4  lb  bird 
»  410  KTS  represents  penetration  by  4  lb  bird 
*  A  of  20  KTS  from  severe  cracking  to  penetration. 

Further,  additional  areas  were  assessed  as  follows: 

From  Fig.  15 

320  -  210  -  110  KTS  (A) 

320  -  230  =»  90  KTS  (A). 

For  F-15  same  areas, 

410  -  110  -  300  KTS 
410  -  90  -  320  KTS. 

Results  summarized  in  figure  10. 

For  the  F-15  canopy  only  one  point  was  tested,  near  the  canopy  arch.  At  160 
KTS,  no  penetration  occurred,  and  at  182  KTS,  penetration  did  occur.  The 
transition  zone  was  estimated  from  180  to  450  using  T-38  percentage  areas  from 
Reference  2.  The  T-38  canopy  value  of  125  KTS  was  represented  for  the  F-15  by 
a  range  from  180  -  450  KTS,  since  the  F-15  canopy  is  thicker  than  the  T-38  and 
it  was  not  felt  that  decreasing  capability  exists  for  the  F-15  as  it  doe3  for 

the  T-38.  The  zone  for  the  T-38  that  is  represented  by  125  KTS  could  have 

easily  been  a  bad  shot.  Therefore,  the  F-15  transition  zones  were  estimated 
as  shown  in  Figure  16.  The  plotted  values  were  bird  penetration  values  from 
Figure  10.  Similar  methods  were  used  to  extrapolate  increased  windshield  and 
canopy  capabilities,  and  these  are  summarized  In  Figures  11  -  14.  Increased  i 

capability  component  strength  distributions  are  summarized  In  Figures  17  and  \ 

18. 

The  above  approach  was  recommended  by  Mr  Blaine  West  from  University  of  * 

Dayton  Research  Institute,  since  it  was  representative  of  realistic  component 
capability.  The  above  component  strength  capabilities  were  estimated  using  * 

the  above  estimation  process.  *  j 


i 

s 

.*> 
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ficuri:  s  profile  OF  F-15F  and  tf  windshield  and  canopy 


300 


CANOPY 


FIG. 18  CRITICAL  AREA  DISTRIBUTION 
NCREASED  CAPABILITY, 500  KT.  NINDSHIElD  AND  350  KT.  CANOPY 


MODEL  APPLICATION 


Since  all  Inputs  have  been  already  derived,  the  number  of  penetrations  on 
a  specific  component,  windshLeld  or  canopy,  can  be  obtained  by  the  product  of: 

//PENETRATIONS  =  OIR  X  FORCE  USAGE  X  AP(V)  X  AP(W)  X  CRITICAL  AREA 

OIR  =»  OPERATIONAL  IMPACT  RATE  PER  106  URR  (15) 

FORCE  USAGE  =  TOTAL  Flight  Hours  for  specific  missLon  air-to-air  or  air — to- 
ground 

A  l’(V)  =  aircraft  velocity  probability 
A  P(W)  =»  bird  weight  probability 
Critical  area  =  Windshield  or  canopy  area  capability 

A  P(V)  and  A  P(W)  values  are  obtained  by  evaluating  the  respective  cumulative 
probability  distribution  junction  for  a  desired  range  (Table  9.0  -  9.6).  The 
critical  area  values  are  obtained  first  by  developing  a  kinetic  energy  matrix 
(Table  10. 0). 

Eqn  14  is  used  to  determine  each  element  of  the  matrix. 

Example : 

For  velocity  range  310  -360  Nominal  335  KTS  and  bird  weight 
4  -  5  Nominal  4.5  LBS 
KE  -  1/2  m  V2 
KE  =  22,344  Ft  Lbs 

For  velocity  range  160  -  210  Nominal  185  RTS  and  bird  weight 
1-2  -*  Nominal  1.5  Lbs 

KE  =  2271  Ft  Lbs 

From  Table  10.0  for  each  kinetic  energy  value  (each  element),  figures  14  -  18 
are  used  to  obtain  proportion  of  component  critical  area.  Table  10.1  and  10.?. 
Lists  the  results  in  the  12  by  15  matrix,  for  the  present  capability 
windshield  and  canopy.  The  use  of  eqn  15  obtains  the  number  of  penetrations 
that  exceed  the  kinetic  energy  of  the  specific  subsystem.  A  summation  is  then 
taken  of  all  incremental  bird  penetration  values  from  the  12  by  13  matrix,  and 
a  total  number  of  penetrations  for  a  specific  usage  profile  is  obtained. 
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AIR/AIR  +  INST/NAV  +  All  Missions  Flown  (except  AIR  TO  GROUND) 

CONUS  &  EUROPE 


CUM.  PROB 

VELOCITY  (KTS) 

VELOCITY  ( KTS) 

AP(V) 

.0001 

116 

110  -  160 

.0948 

.  0949 

160 

160  -  210 

.1501 

.2450 

210 

210  -  260 

.1950 

.4400 

260 

260  -  310 

.2040 

.6440 

310 

310  -  360 

.1660 

.8100 

360 

360  -  410 

.1080 

.9180 

410 

410  -  460 

.0540 

.9720 

460 

460  -  510 

.0210 

.9930 

510 

510  -  560 

.0060 

.9990 

560 

560  -  610 

.0010 

1 . OOOO 

610 

610  -  660 

0 

1 . OOOO 

660 

660  -  710 

0 

1 . OOOO 

710 

F(V) 


-5.793168  +  2. 49687X10” 1  V  -  4.296046X10-3 
+3,774985X10”-’  V3  -  1.792476X10-7  v4 
+4. 396 38 7X10” 10  V5  -  4.371569X10"13  v6 
,  ,  ,  V-40  .2.81  . 

1  “  e^P(-(  307^40 >  > - 


- FOR  116<_V<200 

- FOR  200<  V<  720 


vtot=  279.24  Knots 


TABLE  9.0  F-15  PRESENT  FLEET 
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A I  It  TO  GROUND  MISSION 
CONUS  &  EUROPE 


!UM  PROB 

VELOCITY  (KTS) 

.0002 

110 

.0219 

160 

.0476 

210 

.0587 

260 

.0689 

310 

.1315 

360 

.3220 

410 

.6260 

460 

.8880 

510 

.  9880 

560 

1.0000 

610 

VELOCITY  (KTS)  AP(V) 

110  -  160  .0217 

160  -  210  .0257 

210  -  260  .0111 

260  -  310  .0102 

310  -  360  .0626 

360  -  410  .1905 

410  -  460  .3040 

460  -  410  .2540 

510  -  560  .1000 

560  -  610  .0120 


F(V) 


-4.276016X10-2  +2 
+4.64907 1X1 O'7  VJ 
+3. 569492X1 0_1 2  V5 
,  ,  /V  -  U0_ 

i  expt,-c46r_l  l0 


.515055X10-3  V 
-1.908740X10"9  V* 

-2.406103X10"15 

6.01  . 


■) 


5.1 

v6 


7 6007X10“ 3  V2 


-  FOR  110<V<380 


- FOR  380"  V  610 


vtot*  '*27*90  Knots 


TABLE  9.1  F-15  RAPID  DEPLOYMENT  FORCE 
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AIR  TO  AtR  +  INST/NAV  +  All  Missions  Flown  (except  AIR  TO  GROUND) 

CONUS  &  EUROPE 


CUM  PROB 

VELOCITY  (KTS) 

VELOCITY  (KTS) 

A  P(V) 

.0001 

116 

110  -  160 

.0948 

.  0949 

160 

160  -  210 

.1501 

.2450 

210 

210  -  260 

.1950 

.  4400 

260 

260  -  310 

.2040 

.  6440 

310 

310  -  360 

.1660 

.  8 100 

360 

360  -  410 

.1080 

.9180 

410 

410  -  460 

.0540 

.9720 

460 

460  -  510 

.0210 

.  9930 

510 

510  -  560 

.0060 

.9990 

560 

560  -  610 

.0010 

1 . 0000 

610 

610  -  660 

.0000 

1  . 0000 

660 

660  -  710 

.0000 

1  .  oooo 

710 

-5.793168  +  2.49687X10"1  V 

-  4. 296046X1  O'*3  V2 

F<V)  - 

+3.774985X10-5  v3  _  1.792476X10"'  V'* 
+4.396387X10"10  Vs  -  4. 37 1569X10" 1 3  V"  _J 

- FOR  1 1 6_V;_200 

.  ,  ,  V-40  .2.81  . 

1  exp(-( 307-40'  ' — 

_ 

-  POR  200 '  V  720 

VtoC=  279.24  Knots 


TABLE  9.2  F-15  RAPID  DEPLOYMENT  FORCE 
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AIK 

TO  AIK  9-  T.NST/NAV  +  All  Miss 

tons  Flown  (except 

AIK  TO  GROUND) 

—  :  : 

CONUS 

f,  EUROPE 

UM  PKOB 

VELOCITY  (KTS) 

VELOCITY  (Kl'S) 

AP(  V) 

.0001 

116 

110  -  160 

.0948 

i 

.0949 

160 

160  -  210 

.1501 

.2450 

210 

210  -  260 

.1950 

.4400 

260 

260  -  310 

.2040 

'  >  / 

.  0440 

310 

310  -  360 

.1660 

.8100 

360 

360  -  410 

.1080 

.9180 

410 

410  -  460 

.0540 

.9720 

460 

460  -  5.10 

.0210 

. ; 

.9930 

510 

510  -  560 

.0060 

.9990 

560 

560  -  610 

.0010 

1.0000 

610 

610  -  660 

.0000 

L . 0000 

660 

660  -  710 

.0000 

1.0000 

710 

1 

**  **,l 

'••'.If?  Nit« 

— 

— 

-5.793168  +-  2.49687X10-1  V  -  4. 296046X10- 3  V2 
+3.774985X10~rj  V3  -  1.792476X10-7  V4 

-  FOR  Ufr  V-  200 

"  .  ■  K.  • 

(V)  - 

•F4.396387X10-10  V5  --  4. 371569X10-!3  Vf> 

“  -  :■  ^*1 

.  ,  .  V-40  .2.81  . 

1  ~  exP(“(i(v/_To)  >  . 

Vtot=  279.24  Knots 


TABLE  9.)  F-15  DUAI.  ROLE  FIGHTER 


AIR  TO  GROUND  MISSION 


CONUS  &  Kit  KOI*!' 


!UH  I’ROB 

VELOCITY  (KTS) 

VELOCITY 

(  KTS) 

A  P(  V) 

.0032 

110 

110 

- 

160 

.0431 

.0463 

160 

160 

- 

210 

.0159 

.0522 

210 

210 

- 

260 

.0153 

.0775 

260 

260 

- 

310 

.0195 

.0970 

310 

310 

- 

360 

.0125 

.  1095 

360 

360 

- 

410 

.037  5 

.1470 

410 

410 

- 

460 

.2470 

.3940 

460 

460 

- 

410 

.4390 

.8330 

510 

510 

- 

560 

.1630 

.9960 

560 

560 

- 

610 

.0040 

1.0000 

610 

FOR  IKkVoUO 
FOR  440.-VJ>20 

Vt0[=  445,91  Knots 


TABLE  9.4  F-15  DUAL  ROLF  FIGHTER 


F(V)  = 


+8.284557X10-5  -5.55U1XI0"3  v  +1 . 172276X10"4  V2 
•-9.165418X10-7  V3  +3.493164X10-9  V4 
-6.493033X10-12  V5  +4.716360X10**15  V6 
,  ,  ,V  -  110  .9.5321. 

1  -  e*P(-(wX3-f10)  ' - 


)9.532l) 


CONUS 


TOUT  (LBS) 

A  F(W) 

CUM.  PROB 

WEIGHT 

0  -  1 

.5149 

0 

0 

1  -  2 

.1631 

.5149 

1 

2  -  3 

.1074 

.6780 

2 

3  -  4 

.0668 

.7854 

3 

4-5 

.04.38 

.8522 

4 

5-6 

.0295 

.8960 

5 

6  -  7 

.0204 

.9255 

6 

7  -  8 

.0144 

.9459 

7 

8-9 

.0103 

.9603 

8 

9-10 

.0075 

.9706 

9 

10  -  11 

.0054 

.9781 

10 

11  -  12 

.0040 

.9835 

11 

12  -  13 

.0030 

.9875 

12 

13  -  14 

.0022 

.9905 

13 

14  -  15 

.0017 

.9927 

14 

.9944 

15 

,  „  ,  W  >.4145 

1  "  exp  <-<2.2130' 


FOR  0<W<.9 


V'(W) 


.  1688  W  +  .3461 


-  FOR  ,9<W<  1.7 


,  ,  ,  W  . .7550 

1  exp(-<Y;6955> 


-  FOR  1 .7<W< 1 7 


Wtot»  V.85  Lh. 


TAIU.K  9.5  BIRD  WEIGHT  DISTRIBUTION 


EUROl’li 


rciiT 

( »»B8) 

A  (>(W) 

0 

.. 

1 

.7  793 

i 

.. 

*> 

.1022 

?. 

- 

3 

.0449 

3 

- 

4 

.0243 

4 

- 

5 

.0147 

3 
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.0094 

6 

7 

.0064 

7 

8 

.0045 

8 

- 

9 

.0032 

9 

•o 

.0024 

10 

- 

1  l 

.0018 

1  1 

- 

12 

.0014 

12 

- 

13 

.0010 

13 

- 

14 

.0008 

14 

- 

15 

.0007 

CUM.  l’ROH 

I) 

.7/9:) 

.8815 

.9264 

.9507 

.9654 

.9748 

.9812 

.9857 

.9889 

.9913 

.9931 

.9945 

.9955 

.9963 

.9970 


WHICHT  (L.KS) 
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•> 
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8 
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10 
1  1 
1? 

13 

14 
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TABU;  9.6  BIRD  .JK1CI1T  DISTRIBUTION 
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table  10.2  canopy  critical  area/total  area,  present  capability 


USB  OF  COMPUTER  PROGRAM  FOR  MODEL  APPLICATION 


Appendix  A  lists  an  Interactive  computer  program  which  was  used  to  obtali 
the  number  of  bird  penetrations  for  each  component  varying  tH,  Program  BAAPP 
evaluates  the  present  fleet  air-to-air  mission.  The  Inputs  for  the  program 
are  In  line  number  140  -  740  and  can  he  changed  to  evaluate  the  RIM)  and  IMIF 
a  I  r-to-al r/af r-to-ground  missions.  In  the  above  program  lines,  the  following 
are  Input  parameters: 

(140)  PWC  -  A  P(W)  bird  weight  probability  for  CONUS.  From  Table  9.5  (Woo* 
not  change  when  evaluating  F-15  Present  Fleet,  F-15  RPD,  F-15  l)RF). 

(160)  PWE  -  A  P(W)  bird  weight  probability  for  Europe.  From  Table  9.6  (Doe; 
not  change  when  evaluating  F-15  Present  Fleet,  F-15  RPD,  F-15  ORF). 

(210)  PVI  -  A  P(V)  velocity  probability.  From  Table  9.0  -  9.5  (changes  for 
air-to-air  and  air-to-ground  missions) 

(250  -  740)  Windshield  and  Canopy  Critical  area  Proportions.  From  Table 
10. 1  and  10.2  for  present  Windshield  and  Canopy  Capability. 
(These  values  change  for  increased  capability.  By  using 
Figures  17  and  18,  a  simitar  Table  to  10. 1  and  10.2  can  he 
obta  tned ) . 

Appendix  B  contains  the  results,  after  running  the  program  BAAPF,  for  the 
F-15  present  fleet  (air-to-air  mission).  The  first  matrix,  expected  number  of 
blrdstrikes  (uncorrected),  Is  not  corrected  for  critical  area.  The  second 
matrix,  expected  number  of  blrdstrikes,  is  corrected  for  critical  area,  and 
the  summation  of  all  elements  of  this  matrix  is  the  total  number  of  penetra¬ 
tions  for  a  specific,  component,  theatre,  0IR,  Force  Usage,  or  mission.  The 
Appendix  C  program  lists  only  the  results,  the  total  number  of  bird  penotra- 
t ions . 


EVALUATION  OF  RESULTS 


To  dace,  the  F-15  fleet  has  had  one  bird  penetration  on  the  windshield, 
and  two  Incidents  resulting  In  two  cracked  windshields  and  one  cracked  canopy. 
Of  these  two  incidents,  no  penetration  into  the  cockptt  hns  occurred.  Tahl** 
II  summarizes  the  results  of  this  model's  simulation  of  the  P-15  present  fleet 
baaed  on  the  number  of  hours  flown  to  date  In  Conus  and  Europe  and  on  tin*  pre¬ 
sent  windshield/canopy  capability.  It  Is  Important  to  note  that  any  value 
other  than  an  integer  value  for  a  penetration  is  meaningless;  however,  since 
this  simulation  Is  a  numerical  solution,  fractional  values  will  occur.  Por 
example,  for  Conus,  the  model  predicts  0.45  bird  penetrations  on  the  canopy 
and  to  date  there  have  been  none.  Since  this  is  a  numerical  solution,  one  can 
imagine  a  counter  starting  when  the  P-15  fleet  started  flying.  The  0.45  pene¬ 
trations  implies  we  have  not  had  a  penetration  on  the  canopy  to  this  date,  but 
given  more  time,  we  will. 

The  model  simulated  past  history  very  well,  the  same  methodology  was  used 
to  evaluate  the  F-15  RPD  and  DRF.  The  results  are  presented  in  Figures 
19  -  70  for  different  missions  with  varying  Ts  and  force  usage  (flight  hours). 
A  realistic  usage  profile  is  summarized  In  Table  12.0  -  12.7  for  the  RPI)  and 
DRF  over  a  15  year  flying  period.  These  graphs  can  also  be  used  to  evaluate 
other  and  possibly  more  realistic  usage  profiles. 

Evaluation  of  Table  12.7  Indicates  that,  even  with  a  500  kt  windshield  and 
350  kt  canopy,  bird  penetrations  Into  the  cockpit  will  still  occur;  however, 
there  Is  little  that  can  be  done  to  prevent  an  8  lb  bird  from  penetrating  at 
480  knots.  To  dissipate  the  energy  resulting  from  an  impact  of  such  a  severe 
condition  could  require  structural  considerations  that  would  Interfere  with 
the  capability  and  function  of  the  aircraft,  as  well  as  make  the  design  cost 
prohibitive  based  on  latest  state  of  the  art  technology  In  transparency 
design. 

The  number  of  penetrations  that  are  represented  in  either  the  tables  or 
the  graphs  do  not  imply  aircraft  losses.  Based  on  McAir's  Initial  analyses  of 
historical  data,  one  out  of  every  three  penetrations  resulted  in  an  aircraft 
loss,  fatality,  or  pilot  major  Injury.  This  ratio  was  used  to  assess  the 
results  In  terms  of  aircraft  losses,  fatalities  or  pilot  major  injuries. 
However,  one  should  note  that  this  type  of  ratio  Is  based  on  a  sra.il  I  data 
sample  and  has  not  been  verified  against  all  existing  data  from  which  a  more 
realistic  ratio  might  be  Inferred. 

Figures  19  through  70  can  be  used  to  predict  the  number  of  bird  penetra¬ 
tions  for  a  change  In  Ts  if  the  force  usage  shown  in  Tables  12.0,  12.4,  and 
12.5  remains  the  same.  From  Table  12.3,  it  Is  found  that  the  F-15  RPI)  air  to 
ground  missions  predicts  9,5  windshield  bird  penetatlons  for  a  Ts  of  0.805ft 
and  the  present  windshield  capability  In  CONUS.  What  Is  the  expected  number 
of  windshield  bird  penetrations  If  the  time  spent  In  the  bird  threat  environ¬ 
ment  is  reduced  to  a  Tg  of  0.6?  First  find  the  right  graph.  In  this  case. 
Figure  27,  HF-15  RPD,  Present  Windshield  Capability,  CONUS,  Air  to  Ground.** 
Then  from  Table  12.0,  "F-15  RPD  Force  Usage,”  the  number  of  flight  hours  for 
the  CONUS  and  air  to  ground  case,  141,750  hours.  Is  used  In  Figure  27.  Since 
Figure  27  doesn't  have  a  line  for  a  Tg  of  0.6,  the  value  Is  Interpolated 
between  Ts  of  0.5  and  0.7.  From  the  graph: 


Windshield  Bird 

At  Flight  Hours  Tg  Penetrations  (WBP) 

141,750  .5  WBP. 5  =  5.6 

141,750  .6  WBP. 6  =  ? 

141,750  .7  WBP. 7  =  8.3 

Using  a  ratio  to  interpolate 

WBP. 7  -  WBP. 6  WBP. 7  -  WBP. 5 

Ts.7  ”  ts.6  ts.7  ”  *^s.5 

Solving  for  WBP. 6 
WBP. 6  =  - 

WBP. 6  *>  6.95  windshield  bird  penetrations 

CONCLUSIONS 

A  statistical  simulation  was  used  to  evaluate  the  F-15  present  windshield 
and  canopy  capability  with  respect  to  the  number  of  bird  penetrations  into  the 
cockpit.  The  results  of  this  simulation  as  applied  to  the  present  fleet, 
recreated  past  history  very  precisely.  Similar  methodology  was  used  to  evaluate 
the  F-15  RPD  and  DRF.  The  results  Indicate  that  the  model,  using  the  assump¬ 
tions  and  parameters  of  this  analysis,  predicts  that  the  present  windshield  and 
canopy  need  improvement  for  the  air  -to-ground  mission  of  the  RPD  and  DRF. 


BIRD  PENETRATIONS* 


MODEL  SIMULATION  HISTORICAL 


CONUS 

EUROPE 

CONUS 

EUROPE 

Windshield 

1.10 

0.08 

1 

0 

Canopy 

0.45 

0.62 

0 

0 

Hours  Flown 

473008 

136576 

473008 

116576 

Ts 

.1198 

.2307 

.1198 

.2107 

*with  present  windshield  and  canopy  capability 


TABLE  11  F-15  PRESENT  FLEET  AIR  TO  AIR  MISSION 


Number  of  aircraft  available  »  150  aircraft 

Average  number  of  hours  flown  per  month  »  25  hours 

The  total  number  of  flying  hours  fpr  15  years  Is: 

15  years  X  150  A/C  X  25  hrs  X  12  months  «  675,000  hrs 

month  year 

Aircraft's  distribution:  1  -  Conus  =  70%  of  aircraft 

2  -  Europe  =  30%  of  aircraft 

The  mission  mix  is:  l  -  30%  air  to  ground  (A/G) 

2  -  70%  air  to  air  (A/A) 

Mission  mix  data  obtained  from  11Q  TAC  message  221715%  duly  1982 

For  15  Years  of  F—  1 5  Ri’D  Force  Usage 


Combined 

Conus  (70%) 

Europe  (30%) 

Ts 

Conus 

Ts 

Europe 

Total  hours 

675,000 

472 , 500 

202,500 

A/G  Hours  (30%) 

141,750 

60,750 

.8056 

.8056 

A/A  Hours  (70%) 

330,750 

141,750 

.1198 

.  2  307 

TABLE  12.0  F-15  H PI)  FORCE  USAGE 
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*  "V.  ■v. 


WINDSHIELD  BtRD  PENETRATIONS 


Present 

Capability 

450  Knots 
Capability 

500  Knots 
Capability 

Conus 

.75 

.15 

.09 

Europe 

.08 

.016 

.009 

CANOPY 

BIRD  PENETRATIONS 

Present 

Capability 

300  Knots 
Capability 

350  Knots 
Capability 

Conus 

.30 

.11 

.07 

Europe 

.65 

.19 

.05 

TABLE  12.2  F- 

15  RPD  AIR  TO  AIR  MISSION 

WINDSHIELD  BIRD 

PENETRATIONS 

Present 

450  Knots 

500  Knots 

Capability 

Capability 

Capability 

Conus 

9.50 

3.20 

2.40 

Europe 

.60 

.17 

.12 

CANOPY  BIRD 

PENETRATIONS 

Present 

300  Knots 

350  Knots 

Capability 

Capability 

Capability 

Conus 

2.60 

1.40 

1.05 

Europe 

3.50 

1.50 

1.0 

TABLE  12.3  F-15  RPD  AIR  TO  GROUND  MISSION 
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Production  schedule  obtained  from  HQ  USAF  message  0822007.  Apr  83 


Years  of 
Production 

Number 
of  A/C 

Cumulative 

Number  of  A/C 

Flight  iour9 
per  car 

1st 

16 

16 

4,800 

2nd 

72 

88 

26,400 

3rd 

72 

160 

48,000 

4th 

72 

232 

69,600 

5th 

72 

304 

91,200 

6th 

60 

364 

109,200 

7th 

36 

400 

120,000 

Total  flight 

hours  for  the 

first  7  years: 

447,648 

Sample  Calculation:  for  the  first  year: 

16  A/C  X  25  flight  hours  X  12  month  »  4800  flight  hours 
month  year 

For  the  remaining  8  years,  the  number  of  flight  hours  Is: 

400  A/C  X  25  flight  hours  X  12  month  X  8  years  =  960,000  hours 
month  year 


The  total  flight  hours  for  15  years  is: 

447,648  hours  +  960,000  hours  =*  1,407,648  hours 


TABLE  12.4  F-15  DRF  FORCE  USAGE 


Aircraft  distribution:  1  -  Conus  **  70%  of  aircraft 

2  -  Europe  -  30%  of  aircraft 


The  mission  mix  Is:  1  -  75%  Air  to  ground  (A/G) 

2  -  25%  Air  to  air  <A/A) 


For  15  years  of  F-15  DRF  Force  Usage 


Combined 

Conus 

Europe 

Ts 

Conus 

Ts 

Europe 

Total  Hours 

1,407,648 

985,353 

422,294 

A/G  Hours  (75%) 

739,015 

316,720 

.6360 

.6360 

A/ A  Hours  (25%) 

246,338 

105,573 

.1198 

.2307 

TABLE  12.5  F-15  DRF  FORCE  USAGE 


WINDSHIELD 

BIRD  PENETRATIONS 

Present 

450  Knots 

500  Knots 

Capability 

Capability 

Capability 

Conus 

.55 

.10 

.06 

Europe 

.06 

.01 

0 

CANOPY  BIRD  PENETRATIONS 

Present 

300  Knots 

350  Knots 

Capability 

Capability 

Capability 

Conus 

.AO 

.13 

.08 

Europe 

.80 

.25 

.15 

TABLE  12.6  F-15  DRF  AIR  TO  AIR  MISSION 


WINDSHIELD  BIRD  PENETRATION 


Present 

450  Knots 

500  Knots 

Capability 

Capability 

Capability 

Conus 

45.5 

17.0 

10.9 

Europe 

2.8 

.9 

.6 

CANOPY  BIRD 

PENETRATIONS 

Present 

300  Knots 

350  Knots 

Capability 

Capability 

Capability 

Conus 

19.2 

10.7 

8.9 

Europe 

25.5 

11.5 

8.5 

TABLE  12.7  F-15  DRF,  AIR  TO  GROUND  MISSION 


>000 


urn 


CANOPY  BIRO  PENETRATIONS 
,26  F~15  RPD, PRESENT  CANOPY  CAPABILITY 
EUROPE, AIR  TO  AIR  (0-5000  FT.  AGL) 


WINDSHIELD  BIRD  PENETRATIONS 

-15  RPD  PRESENT  WINDSHIELD  CRPPBILITY 


ship: 

5000 


GOO 


dUO 


AdO 


0  KT.  WINDSHIELD  CRPRB 
GROUND  fQ-BOOO  FF.  RGL 


101 


CRNOPY  BIRD  PENETRATIONS 
G,  38  F-i5  RPD/300  KT.  CANOPY  CAPABILI 
EUROPE,  AIR  TO  GROUND  r0-5000  FT.  AGU 


00b-  0 J  yitj  01  y 
90N1M  0»  009  ‘  Q 


GOh 

NUJ 


.  51  F-15  ORF, PRESENT  WINDSHIELD  CfTPRBILITY 
CONUS, RIR  TO  GROUND  (0-5000  FT.  RGL) 


CRNOPY  BIRD  PENETRATIONS 
-15  DRF, PRESENT  CANOPY  CAPABIL 
AIR  TO  GROUND  (0-5000  FT.  AGl) 


WINDSHIELD  BIRD  PENETRRTIONS 
53  r- 15  DRF,  PRESENT  WINDSHIELD  CflPRBI 
EUROPE, RIR  TO  GROUND  (0-5000  FT.  RGL) 


i  mu'Q 
fr'jyQ  S 


PfiB 

RGL 


CANOPY  BIRD  PENETRATIONS 
-15  DRF, 300  KT.  CANOPY  CAPABIL 
E,  AIR  TO  AIR  (0-5000  FT.  AGL1 


FIG.  59  F- i5  DRF, 450  KT.  WINDSHIELD  CfiPRBILITY 
CONUS, ftIR  TO  GROUND  (0-5000  FT.  RGL) 


Y 


0.0  0.5  i.O  1.5  2. Cl  2.5 

WINDSHIELD  BIRD  PENETR 
FIG.  61  F-15  DRF.450  KT.  WINDSH; 
EUROPE, RIR  TO  GROUND  f0~50C 
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NDSHIELD  BIRD  PENETRATIONS 
:  DRF..500  KT.  WINDSHIELD  CAPABIL 


CANOPY  BIRD  PENETRATIONS 
FIG.  64  F-15  QRF,35G  KT.  CRNOPY  CAPABILITY 
CONUS, AIR -TO  RIR  (0-5000  FT.  RGU 


■5F-15  DRF,500  KT.  WINDSHIELD  CRPRBILITY 
EUROPE,  AIR-  TO  RIR  (0-5000  FT.  RGU 


EUROPE, AIR  TO  GROUND  (0-5000  FT.  RGL) 


i  000000 1  XJ  S^jnOH  1H9HJ 
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5.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0 

CANOPY  BIRD  PENETRATIONS 

FIG. 70  F-15  DRF,350  KT.  CANOPY  CAPABILITY 
EUROPE, AIR  TO  GROUND  IQ-5000  FT.  AGL) 


APPENDIX  A 


100-  PROGRAM  BAAPK { INPUT > OUTPUT) 

110-  DIMENSION  PW<15>,PMC<13),PWE<13>.8WGT<15> 

120-  DIMENSION  VEL < 12) r PV ( 1 2 ) , P01 < 12 ) . ENOB (1 2 . 13 > 

130*  DIMENSION  ACA< 12, 13> . ACAM< 12, 15) . ACAC< 12. 13) 

1 40*  DATA  PWC/. 5149, , 1631 ,. 1074, .0668, .0438, .0295, .0204, ,0144, 

1 50“  1.0103.  ,0073, .0054, .0040, .0030, .0022, .0017/ 

1 60-  DATA  PWE/ . 7703 , . 1 022 , . 0449 , . 0243 ,.0147,. 0094 , . 0064 , . 0043  , 

170-  1.0032, .0024, . 0018 , . 001 4 , . 0010 , .0008, .0007/ 

180-  OAT 4  BHGT/. 3, t  5, 2, 5-. 1.5, 4.3, 5. 3,6 .5, 7. 5, 8. 3, 9. 3, 10.5, 

190-  111.3,12.5,13.3,14.3/ 

200-  DATA  VEL/ "11 0-160" , "160-210" , "210-280" , "280-JlO" , "310-360" , 

210-  1 "360-410" , "410-460" , "460-510" , "310-360" , "360-610" , 

220-  1"G 10-660", "S60-710"/ 

230-  DATA  PV1/.0S48, .  1501 , . 1930,  .2040, . 1660, . 1080, .0340, .0210, .0060, 


240- 

1.0010 

,0.0, 0.0/ 

250-C 

MINDSHIELO  AC/A 

280- 

DATA 

( ACAM ( 1  r  J  > 

,J»1,19>/0.,0.,0.,0.,0.,0.,0.,0., 

270- 

10. ,0. 

i O i iO»  /  0  #  r 

0 .  ,  0 .  / 

280- 

DATA 

( ACAM (  2  ,  J  ) 

» J  3 1 , l3>/0. ,0. ,0. ,0. ,0. ,0. ,0. ,0. , 

280- 

10. ,0. 

, 0. , .238, . 

439, .456, .473/ 

300- 

DATA 

( ACAM ( 3  ,  J ) 

r.hl, 13>/0. ,0. ,0. ,0. ,0. ,0. ,0.0, .432, 

310- 

1.453, 

.407, .314, 

.342,1.0,1.0, 1 ,0/ 

320- 

DATA 

( ACAM ( 4  ,  J ) 

,  J-l , 13)/C. ,0. ,0. ,0. , . 184, .448, .488, .529, 

330- 

11. .1. 

r  i  4  (  1  «  »  1  .  r 

l.,l./ 

340- 

DATA 

( ACAW ( 3  ,  J ) 

,J-l,13)/0. ,0. ,0. , .236, .477, .333, 1 . 0, 1 .0, 

330- 

11  .  ,  1. 

,  i  .  ,  1  .  ,  1 . , 

1. ,  1./ 

360- 

data 

<  ACAW ( 6  ,  J ) 

,  J-l  ,13>/0. ,0, . .! 96,  .484,  .357, 1.0,  1 .0, 1 .0, 

370- 

1  1  .  ,  1  . 

r  ).  •  1 1  n  r  i  «  r 

1.  ,1,/ 

380- 

DATA 

( ACAM  <  7  ,  J ) 

i*  J  *  1 » 13)/0»  rO*  t  •  46 1 1  •  353  *i«0»l«0rl  .Or  1  »0r 

390- 

11. ,1. 

^  1  4  r  1  ,  4  1  ,  ) 

!.,!./ 

400-  DATA  (ACAU(G,.I),  J«i,15)/0.,  0.0,  .518,1. ,  1 . ,  1 . ,  1 , ,  1 . , 

410-  1 1 . , 1 . , i  . , 1 .  ,  1. ,  1.  ,  1 ,/ 

420-  DATA  ( ACAM <9, J) , J-l . 13>/0. , .440. x. 0,1. 0,1. 0,1. 0,1. 0,1.0, 

430—  11, pl.,1.,1. , 1 . , 1 .  ,1./ 

440-  DATA  ( ACAM  (  1  0 ,  J )  .  ,J  - 1  .  1 3 )  /O .  , .  48 1,1. 0,1. 0,1. 0,1. 0,1. 0,1.0, 

430-  ll.rl.rl,,l«,l<,l.,l./ 

460*  DATA  ( ACAW( 11 , J > , J»1 , 13) /O. ,.327,1.0, 1.0, 1.0, 1.0,1. 0,1.0, 

470-  11.  ,1.,  1 1  •  ,1.,1.,1./ 

480-  DATA  <ACAW< 12, J) . J-i , 15J/0.0, 1 . , 1. , l . , 1 . , 1 . , 1 . , l .0, 

490-  U.,1.,1.,1.,1.,1.,1./ 

300-C  CANOPY  AC/A 

310-  DATA  <  ACACU  .  J)  ,-1-1 ,  l3)/0.  ,0.  ,0.  ,0.  ,0.  ,0.  ,0.  ,  .207, 

520-  1.302, .398, . 433, , 3B0 , . 684 , . 779 , .805/ 

330-  DATA  <  ACAC  <  2 . J  > . J  « 1 , 13 ) /O . . 0 . , 0 . , 0 . , . 298 , . 477 , . 636 , . 802 , 

340-  1  .814,  .626, . 833 , . 631 » . 063 , . 873 . . 888/ 

530-  DATA  ( ACAC ( 3 , J ) .  1-1 , 13>/0. ,0. , .214, .303, .792, .819, .838, .83B, 

560-  1.878,  .898, .318, .937, .937, .977, .997/ 

370-  DATA  <  ACAC (  4 .  J )  ,  J«f. ,  19)/0.  ,0.  ,  .334,  .812,  .841 ,  .870,  .899,  .928, 

380-  1.937,  .906,1.0,1.0,1.0,1.0,1.0/ 

390-  DATA  < ACAC  <  5 , J  > , J«i . 131/0. , .373, . 810 . . 83 1 , . 391 , . 93 1 , . 97 1 , 1 . 0 , 

600-  11.0,1.0,1.0,1.0,1.0,1.0,1.0/ 

610-  DATA  < ACAC<6. J)  ,  J«1 , 15>/0.  ,  .633, .843, .096. .949, 1.0, 1.0, 1.0, 

620-  11.0,1.0,1.0.1.0,1.0,1,00,1.0/ 

630-  DATA  <  ACAC ( 7, J ) , J-l , 13) /O. . .812, .879, .947, 1.0, 1.0, 1.0, 1.0, 

640-  11.0,1.0,1.0,1.0,1.0,1.0,1.0/ 

650-  DATA  < ACAC(9. J) . J»1 , 13)/0. . .836, ,921 , 1 .0, 1 .0, l.O, 1 .0. 1 .0, 

660-  11.0,1.0,1.0,1.0,1.0,1.0,1.0/ 

670-  DATA  ( ACAC <9, J) , J- 1 , 15 ) / .241 , .864, .966, 1.0, 1.0, 1.0, 1.0, 1.0, 

680-  11.0,1.0,1.0,1.0,1.0,1.0,1.0/ 

680-  DATA  <  ACAC i 10,3) , J-l , 13) /  .38/, .894, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 

700-  11.0.1.0,1.0-1.0,1.0,1.0,1.0/ 


13/1 


710- 
720- 
730- 
740- 
730- 
780- 
770- 
7G0- 
790- 
800- 
810- 
820- 
830- 
840- 
850- 
860- 
870- 
880- 
890- 
800- 
910- 
920- 
930- 
940- 
930- 
960- 
970* 
980- 
990- 
1000- 
101O* 
1020- 
1030- 
1040- 
1030- 
1060- 
1070- 
1080- 
1090- 
1100* 
1110" 
1120- 
1 1  so¬ 
il  40- 

I  ISO- 
1160- 
1170- 

I I  So¬ 
il  90- 
1200- 
1210- 
1220- 
1230- 
1240- 
1250- 
1260- 
1270- 
1280- 
1290- 
1300- 


DATA  < ACAC (11,J),J- 1,13)/. 347.. 927, 1.0, 1.0, 1.0, 1.0, l, 0,1.0 

11.0,1.0,1.0,1.0,1.0,1.0,1,0/ 

DATA  (ACAC < 12, J) , J- 1,15)/. 718, .962,1.0,1.0,1.0,1.0,1.0,1.0 
11.0,1.0,1.0,1.0,1.0,1.0,1.0/ 
t  PRINT  * , "WINDSHIELD  OR  CANOPY  ?  * 

READ  1010, PART 

IF<  PART .EG. “WINDSHIELD* ) GOTO  810 
IF  <  PART . EQ . "CANOPY* ) GOTO  S30 

PRINT  4, "PLEASE  ENTER  'WINDSHIELD '  OR  'CANOPY'* 

GOTO  t. 

810  DO  820  1*1,12 
DO  820  J-1,15 
820  ACA  <  I ,  J  > -ACAW< I , J  > 

GOTO  5 

930  DO  840  1-1,12 
DO  940  J-1,1S 
840  ACA  < I , J ) -ACAC ( I »  J  > 

3  PRINT  *, -OPERATIONAL  IMPACT  RATE  (  PER  10**6  >  • 

READ  *,OIR 

10  PRINT  *, -CONUS  OR  EUROPE?  -  - 

READ  10 10, WHERE 

PRINT  *," ENTER  FORCE  USAGE  -  “ 

READ  *,FU 

IF < WHERE. EQ.-C0NU8-) GOTO  20 

IF ( WHERE. EQ. -EUROPE") GOTO  40 

PRINT  *. "PLEASE  ENTER  ' CONUS '  OR  '  EUROPE'" 

GOTO  10 

20  DO  30  1-1,13 
30  PW< I ) -PWC< I > 

GOTO  90 

40  DO  30  1-1,13 
30  PH( I > -PWE ( I ) 

60  PRINT  *. -ENTER  X  F-13  •* 

READ  *, P13 

IF < PI 3. EG. 100. )QOTO  70 

PRINT  *, -WARNING  THE  NIX  DID  NOT  TOTAL  100X" 

70  F13-P13/100* 

DO  80  1-1,12 
80  PM(  I ) -F13*PM1 < I ) 

1P13-P13 

PRINT  1 05  0, WHERE, IP13, PART, OiR.FU 
PRINT  10'„3,BWQT 
DO  1-1,12 
DO  SO  J-1,13 

ENOB < I , J  > -0 1 R»FU*PM < I >  *PW  <  J ) 

90  CONTINUE 

100  PRINT  1030 , MEL ( I ) , <  ENOB( I, J), 3- 1,13) 

PRINT  1040 
PRINT  1023, BWGT 
6-0 

DO  120  1-1,12 
DO  110  J-1,13 

ENOB ( I , J ) -ENOB  < I , J  >  *AC A ( I , J  > 

B-EN08 ( I , J )  +8 
110  CONTINUE 

120  PRINT  1030, MEL < I ) , <ENOB( I , J ) , 3-1 , 13) 

PRINT  1060,8 
1010  FORMAT ( A10) 

1020  FORMAT ( 1H1 , 

1  10X.A10, 13,"  F-13  MIX  C".A10,">", 


1310-  2  /  /  •  10X ,  "OIR  -".F16.2,"  /  10##8  FU  =  "»F10.4,"  *  10**6", 

1320-  3  /  /  /  / , 25X t "EXPECTED  NUMBER  OF  BIRDSTRIKE  (UNCORRECTED) - > 

1330-  1025  FORMAT* /,3X. "VELOCITY  " ,25X# "BIRD  WEIGHT  (LBS)", 

1340.  1  /  ,  3X» "  < KNOTS,  %15F8.3> 

1350-  1030  FORMAT ( /  , 3X» AlO, 15F8.4) 

1380-  1040  FORMAT* ////,2SX, "EXPECTED  NUMBER  OF  BIRDSTRIKES  (CORRECTED)") 

1370-  1060  FORMAT  (  /  *  25X  t-  *  THE  TOTAL  EXPECTED  NUMBER  OF  PENETRATI0N8  -  "-F10.4) 

1380-  END 

1380«*E0R 
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APPENDIX  B 


VELOCITY  BIRD  HEIGHT  <LBS> 

(KNOTS)  .300  1.500  2.500  3.500  4.300  5.300  6.300  7.500  8.300  9.300  10.300  11.500  12.300  13.300  14.300 
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EXPECTED  NUMBER  OF  BIRDSTRIXES  (CORRECTED) 


430006  CM  STORAGE  USED 

.906  CP  SECONDS  COMPILATION  TIME 
UIND SHIELD  OR  CANOPY  7CAN0PV 
OPERATIONAL  IMPACT  RATE  (  PER  10*«8  >  -31.23 


0000 


VELOCITV  BIRD  HEIGHT  (LBS) 

(KNOTS)  .500  1.500  2.500  3.500  4.500  3.500  6.500  7.300  8.300  9.300  10.300  1 1  .T500  12.500  13.500  14.500 
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APPENDIX  C 


100*  PROGRAM  RAAPF< INPUT, OUTPUT) 

110*  DIMENSION  PW < 1 3 ) , PWC ( 1 5 >  *  PME  < 1 3 > , BWGT < 1 5 ) 

120*  DIMENSION  UEL< 12  > , PO < 12 ) « PU1 <  12) , ENOB( 12,15) 

130*  DIMENSION  ACA< 12 , 15> , ACAH< 12 , 15 > , ACAC< 12 , l 3 > 

140“  DATA  PWC/. 3149, .1631, .1074,  .0668, .0438, .0293, .0204, .0144, 

130*  1.0103, .0073, .0034, .0040, .0030, .0022,  .0017/ 

160*  DATA  PME/. 7793 ». 1022, .0449, .0243, .0147, .0094, .0064, .0043, 

170*  1.0032, .0024, .0019, .00  1  4, .0010, .0008,  .0007/ 

180*  DATA  BWGT/. 3, 1.3, 2. 3, 3. 3, 4. 3, 5. 3, 6. 3, 7. 3, 8. 3, 9. 3, 10. 3, 

190*  111.3,12.3,13.3,14.3/ 

200*  DATA  V£L/“ 110-160", “160-210“, “210-260", "260-310", H310-3G0M, 

210*  1 "360-410“ , “410-460" , "460-310“ , "310-360" , "360-610" , 

220*  l "610-680", "660-710"/ 

230*  DATA  PV1/.0948, . 1501 , . 1950,  .2040, . 1660, . 1080, .0340, .0210, .0060, 

240*  1.0010,0.0,0.0/ 

2S0-C  WINDSHIELD  AC/A 

260*  DATA  (ACAW( 1 , J) , J-l . 13)/0.  ,0. ,0. ,0. ,0. ,0. ,0. ,0. , 

270*  10. , 0 • , 0 • , 0 . , 0. , 0 • , 0. / 

260-  DATA  ( ACAW<  2 ,J) » J-l » 13) /O.  ,0. ,0. ,0. ,Q. ,0. ,0.  ,0. , 

290*  10. ,0. ,0. , . 258 , . 438 , . 436 , . 473 / 

300*  DATA  ( ACAM<  3, J  > ,J*I,13)/0. , 0 . , 0. ,0. ,0. ,O.,0.O, .432 , 

310*  1.439, .487, .314, .342, 1 .0,1.0, 1.0/ 

320*  DATA  (ACAW<4, J) . J*l. 13)/0. ,0. ,0. ,0. ,.184, .448, .488, .329, 


330- 

11. 

r  1  • 

340- 

DATA 

(ACAM<S,J> , J-l, 131/0. ,0.,0. , .236, .477, .333,1.0, 1 

.0, 

350- 

11. 

r  1  • 

360- 

DATA 

(ACAH<6, J) , J-l , 19) /O.  ,0. , .198. .484, .337, 1.0, 1 .0, 

1.0, 

370- 

11. 

9  1  • 

380- 

DATA 

(ACAN(7,J) . J-l, 13)/0.,0.,. 461,. 353. 1.0, 1.0, 1.0,1 

•  Or 

390- 

11. 

9  1  • 

400- 

DATA 

(ACAH<8,  J) ,  J-l,13)/0.  ,0,0,  .319,1.  ,.l .  ,1. ,  l.,l.  , 

410- 

11. ,1. 

,1. ,1.,1.,1.,1./ 

420- 

DATA 

(ACAM(9. J) ,J-l,13)/0., .440,1.0,1.0,1.0,1.0,1.0,1 

•  Or 

430- 

11. 

»!• 

440- 

DATA 

<  AC  AW  (10,J),.J«l,13)/0.  ,  .481,1.0,1.0,1.0,1.0,1.0, 

1.0, 

450- 

11.  ,1. 

480- 

DATA 

(ACAUUl,  J),  J-l,  13) /O.  ,  .327,  1.0, 1.0, 1.0,  1.0 ,1.0, 

1.0, 

470- 

11. 

«  1 « 

480- 

DATA 

( ACAM < 12, J) . J-l » 13)/0.0, 1. , 1  .  ,1. , 1 .  ,1 . , 1 .  ,1.0, 

490- 

11. 

9  1  • 

500*C  CANOPY  AC/A 

510-  DATA  <ACAC< 1, J) , J-l , 13>/0. ,0. ,0. ,0. ,0. ,0. ,0. , .207, 

520-  1.302, .398, .493, .368. .684, .779, .803/ 

330-  DATA  (ACAC(2, J> , J-l . 13)/0. ,0. ,0. ,0. , .298, .477, .658, .802, 

540-  1.814; .828, .838, .851 , .883, .873. .888/ 

330-  DATA  (ACAC(3, J) , J-l . 13>/0. ,0. » .214, .503, .792, .613, .833, .638, 

360-  1.878, .898, .918, .937, .937, .977, .997/ 

370-  DATA  <ACAC<4, J> , J-l. 13>/0.  ,0. ,.334, .812, .841 , .870, .899, .928, 

360-  1.937, .986, 1.0, 1 .0, 1 . O » 1 .0, 1 . 0/ 

590-  DATA  <ACAC(3, J) , J-l, 13)/0. , .373, .810, .831 , .891 , .931 , .971 , 1 .0, 

600-  11.0,1.0, 1 .0,1. 0,1. 0,1 .0,1.0/ 

610-  DATA  <ACAC(8, J) , J-l , 13)/0. , .653, .843, .898, .949, 1 .0 ,1 .0, 1 .0, 

620*  11.0,1.0, 1 .0,1. 0,1.0, 1 .00,1 .0/ 

630-  DATA  <ACAC(7,J> , J-l , 13>/0. , .812, .879, .947,1.0,1.0, 1 .0, 1 .0, 

640-  11.0,1.0,1.0,1.0,1.0,1.0,1.0/ 

630-  DATA  (ACAC(8, J) , J-l , 15)/0.  , .836, .921, 1 .0,1 .0, 1.0, 1 .0, 1 .0, 

660-  11. 0,1. 0,1. 0,1. 0,1. 0,1. 0,1.0/ 

670-  DATA  <ACAC<9. J) , J-1, 13)/. 241 , .664, .966, 1 .0,1. 0,1. 0,1. 0,1.0, 

680-  11.0,1.0,1.0,1.0,1.0,1.0,1.0/ 

690-  DATA  <ACAC( 10, J) , J-l , 1 3) / . 387, . 894 , 1 . 0 , 1.0, 1.0, 1.0, 1.0, 1.0, 

700*  11. 0,1.0,1.0,1.0,1.0,x. 0.1.0/ 
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710*  DATA  (ACAC!11,J>,J»1. 15)/. 547,. 927, 1.0,1. 0,1. 0,1. 0,1. 0,1.0, 

720*  11.0,1.0,1.0,1.0,1.0,1.0,1.0/ 

730*  DATA  (ACAC!12, J)»J>  1,15)/.  719,.  962, 1.0, 1.0, 1.0, 1.0, 1.0 ,1.0, 

740*  11.0,1.0,1.0,1.0,1.0,1.0,1.0/ 

730-  1  PRINT  *,  "WINDSHIELD  OR  CANOPY  <M  OR  C)  ?  -" 

760-  READ  1010,  PART 

770*  IF!  PART. EO. "W" ) GOTO  810 

780*  IF! PART. EQ. "C" ) GOTO  830 

790*  PRINT  *, -PLEASE  ENTER  'W'  FOR  WINDSHIELD  OR  'C'  FOR  CANOPY ' ■ 

800*  GOTO  1 

810-  810  DO  820  I  >1,12 

820-  DO  820  J> 1 , 15 

830*  820  ACA< I , J) -ACAWC I ,  J> 

840-  GOTO  5 

850*  830  DO  840  I >1,12 

860*  DO  840  J*  1 , 15 

870-  840  ACA ( I , J ) • ACAC ( I , J ) 

880*  5  PRINT  *, "OPERATIONAL  IMPACT  RATE  (  PER  10**6  )  -  • 

890*  READ  *,OIR 

800*  10  PRINT  *, -CONUS  OR  EUROPE?  !'C'  OR  '£') - • 

910*  READ  1010 , NHERE 

320-  IF<  WHERE.  EQ.  HCn  .  OR  .WHERE  •  EQ .  "E" )  00  TO  15 

930-  PRINT  *, -PLEASE  ENTER  'C'  FOR  CONUS  OR  'E '  FOR  EUROPE - - 

940-  GOTO  10 

930*  15  PRINT  *,  -ENTER  FORCE  USAGE - - 

960-  READ  *,FU 

970-  IF ( WHERE. EQ. "E" )  GO  TO  40 

980-  20  DO  30  1-1,13 

990-  30  PH ( I ) -PNC ( I ) 

lOOO-  GO  TO  60 

1010-  40  DO  30  1*1.13 

1020-  50  PW< I ) -PWE < I ) 

1030-  60  PI 5-100 

1040*  IF  ( P15.EQ .100. )GOTO  70 

1050-  PRINT  *,  "WARNING  THE  MIX  DID  NOT  TOTAL  100X" 

1060*  70  F15-P15/100. 

1070-  DO  80  1-1.12 

1080-  80  PV! I ) -F13*PV1 !  I ' 

1090-  IP13>P15 

1100-  DO  100  I • 1 , 12 

1 1 10-  DO  90  J-l  .  13 

1 1 20-  ENOB ( I ,  J ) -GIR*FU*PV  < I >  *PW ( J ) 

1130-  90  CONTINUE 

1140*  100  CONTINUE 

1130-  8-0 

1160-  DO  120  1*1,12 

1170-  DO  110  J*  1 , 13 

1 1 80-  ENOB U , J ) -ENOB < I, J) * ACA ( I ,  J) 

1190-  B>ENOB< I , J )  +8 

1200-  110  CONTINUE 

1210-  120  CONTINUE 

1220*  PRINT  1060, B 

1230-  1010  FORMAT! A10) 

1240*  1020  FORMAT!  1H1. 

1230-  1  10X , A10 , 13 , "  F-15  MIX  (",A10.">-, 

1260-  2  // , 10X, "OIR  *",F10.2»"  /  10**8  FU  -  ",F10.4,"  *  10**6", 

1270-  3  ////,25X, "EXPECTED  NUMBER  OF  BIRDSTRIKE  (UNCORRECTED)") 

1280-  1025  FORMAT </,3X," VELOCITY  -,23X,"BIRD  WEIGHT  !LBS>", 

1290-  1  / , 3X, • ! KNOTS)  ",13F8.3> 

1300-  1030  FORMAT!/, 3X,A10,13F8. 4) 

1310-  1040  FORMAT! ////,23X, "EXPECTED  NUMBER  OF  BIRDSTRIKES  ! CORRECTED )“ ) 

1320*  1080  FORMAT!/ . 23X, "THE  TOTAL  EXPECTED  NUMBER  OF  PENETRATIONS  -  " ,F10.4? 
1330-  END 
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